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ABSTRACT 

The standard cooling flow model has predicted a large amount of a cool gas in the clusters of galaxies. The failure of the Chandra and 
XXM-Newton telescopes to detect a cooling gas (below 1-2 keV) in the clusters of galaxies has suggested that some heating process 
must work to suppress the cooling. The most likely heating source is the heating by AGNs. There are many heating mechanisms, but 
we will adopt the effervescent heating model which is a result of the interaction of the bubbles inflated by AGN with the intra-cluster 
medium(ICM). 

Using the FLASH code, we have carried out ID- time dependent simulations to investigate the effect of the heating on the suppression 
of the cooling in cooling flow clusters. We have found that the effervescent heating model can not balance the radiative cooling and it 
is an artificial model. Furthermore, the effervescent heating is a function of the ICM pressure gradient but the cooling is proportional 
to the gas density square and square root of the gas temperature. 
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1. Introduction 

According to the steady flow assumption of the standard cooling 
flow model, we must find a cool gas and a multi-phase medium 
within the cooling core in clusters of galaxies which are not ob- 
served in any wavebands (X-ray and non X-ray). This is known 
as the cooling flow problem in clusters of galaxies. In others 
words, there is a discrepancy between standard cooling flow 
model and observations (X-ray and non X-ray). This strong dis- 
crepancy is interpreted as either the gas is being prevented from 
cooling by some heating process, or it cools without any spec- 
troscopic signature (Fabian et al. 2001) which is difficult. But, 
we see that this discrepancy between the standard cooling flow 
model and X-ray observations indicates that either the cooling 
in the center of cooling flow clusters must be suppressed by any 
heating mechanism, or the steady flow assumption of the stan- 
dard cooling flow model is not appropriate. Somewhere else, we 
have concentrated in the second point which it is found that the 
steady flow with cooling is impossible, i.e the cooling flow prob- 
lem is due to the wrong steady flow assumption. In this work we 
will concentrate in the first point, heating the cooling flow. 
The failure of the multi-phase model has revived the idea of a 
heating mechanism which can suppress the cooling. There are 
five main conditions for the heating (see Gardini & Ricker 2004, 
for review) : 

1. The heating must be fine tuned and distributed to get the 
smooth observed temperature profile. Too much heating 
would result a outflow from the center region. Too little heat- 
ing is not sufficient to suppress the cooling. Moreover, the 
heating process must be self regulated: the mass flow rate 
triggers the heating and the heating reduces the mass flow 



rate (Bohringer et al. 2002). That mechanism is called a heat- 
ing with a feedback. 

The heating mechanism by AGN must be sporadic because 
the radio activities are not observed in every cooling flow 
clusters. 

The kinetic energy injection must be subsonic and the ICM 
must not be shocked, as observed in general case. The shock 
could compress the gas producing the catastrophic cooling 
much faster. 

The heating mechanism must preserve the observed entropy 
profiles of the ICM, which decrease toward the center. 
The heating must not destroy the metallicity profiles which 
are peaked toward the centers (Tamura et al. 2002; De Grandi 
& Molendi 2001; Irwin & Bregman 2001). 



A giant elliptical or cD galaxy sits at the potential well of ev- 
ery cooling flow cluster (Mathews et al. 2003; Eilek 2004). The 
popular heating mechanisms are a heating by AGN (Binney & 
Tabor 1995; Tabor & Binney 1993; Bruggen & Kaiser 2002), 
thermal conduction, cosmic rays, galaxies motions, magnetic 
field reconnection (Soker & Sarazin 1990) and turbulent mix- 
ing (Kim & Narayan 2003). About 71 percent of the cDs in the 
cooling flow clusters are radio loud compared to only 23 percent 
of non-cooling flow clusters (Burns 1990). This result suggests 
that there a relationship between the AGN activities and the pres- 
ence of cooling flow. 

In some of cooling flow clusters, the recent X-ray observa- 
tions reveal holes in the X-rays surface brightness coincident 
with the radio lobes, known as X-ray cavities or bubbles ( see 
Fabian et al. 2006, as example). The bubbles or cavities are not 
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a universal phenomenon, suggesting a duty cycle. Radio sources 
are not even detected in some cooling-flow clusters. 

The observed metallicity gradients make a constraint on 
the abiUty of baubles to mix the ICM (Bohringer et al. 2004). 
Brighenti & Mathews (2002) have run many ID simulations of 
clusters with various levels of heating. They concluded that the 
best fit of temperature profiles with the real clusters is without 
heating. The AGNs are assumed to inject buoyant bubbles into 
the ICM, which heat the ambient medium by doing work PdV as 
they rise and expand. This mechanism is called the effervescent 
heating model or mechanism (Begelman 2001). In this work, we 
will concentrate on this mechanism, using time dependent hy- 
drodynamics simulations, FLASH code (Fryxell et al. 2000). 



2. Effervescent heating model 

The center AGN is assumed to inflate bubbles of relativistic 
plasma in the ICM. These bubbles will expand as they rise, doing 
PdV work on their surroundings. Assuming a steady state and 
spherical symmetry, the energy available for heating the ICM 
known as the effervescent heating (Begelman 2001; Ruszkowski 
& Begelman 2002) is: 

E oc Phirfyx-^^l^" , 

where Ph is the partial pressure of buoyant fluid inside the bub- 
bles at radius r, and jb is the adiabatic index of buoyant fluid. 
Begelman (2001) and Ruszkowski & Begelman (2002) have as- 
sumed that the pressure of a buoyant fluid is scaled with the ther- 
mal pressure of the ICM; i.e. the ratio between the pressure of a 
buoyant fluid and the thermal pressure of the ICM is a constant 
with radius. In this case: 

E oc Pijiyo-^^l^t, 

In that model, we should note that the heating profile by AGN is 
proportion to the ICM properties, not the AGN itself. The heat- 
ing rate per unit volume is given by: 



-h{r) V • 



AUr^ = -'^'^ [p 



r dlnr 



(1) 



where 'Hv is heating rate (erg cm^^ 5ec."'), P,, is the center 
pressure and h{r) is a normalization function multiplied by a cut- 
off function (1 - e'''^''") as proposed by Ruszkowski & Begelman 
(2002) : 



hir) = 



4nr^ 



q-^ X (1 - e-'''") , 



where is a small cutoff radius (10-25 kpc) and L is the 
Luminosity of the center source which is given by: 



r 

■Jo 



Upon reflection, the function (1 - e^''^''") is not a cutoff function 
as proposed, but it is a shallowness function added to let the 
heating profile shallower more than the original one, working as 
a artificial fine tune. The function e^''^''" is working as a cutoff 



function. As increases, we get a shallower profile. 
The factor q is given by 



1 dlnP 



r d]nr 



(1 - e-'''") dr 



We found that without the minus sing (not in the original paper) 
heating will be negative. This factor ^ is a normaUzation factor 
to let the total heating inside r^j^ be equal to L. By doing this 
normalization, we find that the presence of the term Pg (in equa- 
tion 1) is not essential. 

Finally, the luminosity of the central source L is given by a feed- 
back from mass flow rate (mass accretion rate) 

L = eMfiow-min C 

where e is the accretion efficiency and c is the speed of light. 
Mfiow-min which is given by: 



Mfla 



= (-4^'-mmPV) 



is the mass flow rate in the inner radius r„„„ depending on the 
resolution of simulation. For high resolution simulation, we are 
aware that the mass flow rate is a small quantity which is not 
enough to fuel the AGN but that model is most accurate one 
we know. The actual heating mechanism is still unclear and 
uncertain. 



At the end, we get the heating profile as: 
L „^ dP 



P^ ^(l- e-'l'") 



An qr^ dr 

and the internal energy equation becomes: 

c*//v dP 
^.V(//.v)-- = -...^. 



(2) 



(3) 



where e^, is the X-ray emissivity and is the gas enthalpy per 
unit volume (U^ + P), and is given by: 



5kT, 



2 fi nip 
where pg is the gas density. 

2.1. Problem with the effervescent heating model 

The effervescent heating profile depends on the gradient of ICM 
pressure rather than temperature and it is much steeper than the 
X-ray emissivity. This model needs a fine tuning because the 
velocity in the center must be a negative quantity, otherwise the 
gas could be accumulated somewhere near to the center. The 
gas will flow from outwards and inwards into that region. In 
other words, the feedback can not ensure the fine tuning and can 
cause the velocity in the center to be a positive quantity and in 
the other part of the cooling core to be negative. 

Ruszkowski & Begelman (2002) claimed that their model 
does not need a fine tuned heating. We argue that their model is 
fine tuned by choosing some value of to control the steepness 
of the heating profile, and choosing e accretion efficiency to 
control the amphtude of the heating function. The accretion 
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efficiency is not a constant for any simulation but is allowed to 
change to guarantee an artificial fine tuned heating. 
The scaling of the buoyant fluid pressure inside the bubbles with 
the ICM pressure is a difficult assumption, during the motion of 
bubbles in ICM. The bubbles must expand in some way to keep 
this ratio a constant. The free parameters (rg, r„jax, e, and r,„,„) 
in the effervescent heating model are responsible for the model 
to be artificial and do not reflect any physics for the heating. 

3. Initial conditions 

The total cluster mass or the gravitation acceleration can be de- 
termined from X-ray observations by assuming that the gas is in 
a hydrostatic equilibrium and a spherical symmetry 



g ^ — 



VP 

Pg 



(4) 



where P and pg are the observed gas pressure and density. The 
inflowing gas has nearly no effect on the total cluster mass, then 
we can assume that the value of g nearly does not change since 
the cluster started to cool. 

In cooling flow clusters, the temperature and gas density 
profiles change only in the cooling core region due to cooling, 
but they do not change in the outer part of the cooling core. The 
cooling has no effect on the total clusters masses. We determine 
the initial temperature profile by fitting the observed X-ray 
temperature profile of the outer part of the cooling core with the 
universal temperature profile of Loken et al. (2002). 



T(r) 



= To{l + fl.,- r/r,,,-,.) 



(5) 



where To, a^, and 6 are free parameters. For most clusters, these 
fitting parameters are found to be a ^ - 1,5 , 6 = 1.67, and 
To = 1.31 . Then the gas density is given by 



Pg('") = Pgo X 
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Fig. 1. Initial conditions: the dots are the observed X- 
and gas density profiles (Vikhlinin et al. 2005). The 
initial fitted temperature profile (Loken et al. 2002). 
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Fig. 2. Model A, pure cooling: the simulated gas density and temper- 
ature profiles at different times up to 7.2 G yr. The big dots represent 
the observations. The simulation agrees very will the X-ray observa- 
tion (Vikhlinin et al. 2005). 



where pg is the gas density, pgo and Tg are the gas density and 
temperature at the center. The g{r) is obtained by current X-ray 
observation assuming a hydrostatic equilibrium as in equation 4. 

We will take the observed X-ray gas and temperature profiles 
of the Chandra observation (Vikhlinin et al. 2005) as an initial 
conditions for the cluster A 1991 as in figure 1. 

3. 1. Models 

We have carried out four model simulations assuming the effer- 
vescent heating model. Model A is carried out only with cooling. 
Model B is a model with an effervescent heating and cooling 
and model C is the same as model B but with higher resolution 
(smaller value of r„„„). Finally, the model D has a higher cutoff 
radius rg to get a smoother heating profile, see table 1 . 



4. Results and discussion 

In this section, we describe our four simulations and the results 
will be given as follow: 

4.1. Model A 

The cluster AI99I is observed by the Chandra X-ray observa- 
tory (Vikhlinin et al. 2005). Model A is a time dependent pure 
cooling simulation, there is no heating. This simulation have run 
until time 7.2 Gy where the simulated temperatures and densi- 
ties agree very well with the X-ray observations (Vikhlinin et al. 
2005) (see figure 2). After time 7.2 Gy, the catastrophic cooling 
accrued only within a few cells due to the resolution, reflecting 
boundary conditions and steep potential profile. The heating can 
not help because it will heat the whole cooling core. 
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Table 1. The models 



Model 


Note 


cutoff radius 


max 


e 


Resolution -r„„„ 


Model A 


Pure Cooling 








520 pc 


Model B 


Cooling and Heating 


22Kpc 


2000 Kpc 


.01 


520 pc 


Model C 


Cooling and Heating 


22Kpc 


2000 Kpc 


.01 


78 pc 


Model D 


Cooling and Heating 


40Kpc 


2000 Kpc 


.001 


520 pc 




4.2. Model B 

The mass accretion rate (to fuel the black hole) is taken to be 
equal to the mass inflow rate at the inner radius r„„„. In this 
model, we have set the cutoff radius r„ and heating efficiency 
e at 22 kpc and .01 respectively, see table 1. The resolution in 
the center is about 520 pc which is large than the black hole 
accretion radius. Ruszkowski & Begelman (2002) have chosen 
r,„i„ equal to 1 Kpc to get large mass flow rate which can fuel 
the black hole in the cluster center. We are certain that at small 
radius the mass flow rate is not sufficient to fuel the central black 
hole to stop the cooling. However, that is the only way to study 
the effect of the heating in clusters properties. The heating pro- 
file is much steeper than the X-ray emissivity which makes it 
difficult for the heating profile to balance the radiative cooling. 
In figure 3, we plot the heating rate and cooling rate ratio. As in 
figure 4, the cluster center is overheated and there is no balance 
between the heating and cooling; that is clear in the temperature 
profile. Moreover, inside 10 kpc, the gas density peaks near to 
center not in the center which is not observed. The gas becomes 
a clumpy structure. There is a outflow from the cluster center 

4.3. Model C 

In model B, the inner radius (resolution) is about 520 pc which is 
larger than the black hole accretion radius. This model is similar 
to model B but the resolution is higher to get smaller accretion 
radius (inner radius). We have set the resolution or inner radius 
equal to 78 pc which can not be achieved by three dimensions 
simulations. The result is that the heating is not enough to bal- 
ance the radiative cooling because the mass flow rate at a small 



radius is very small . Of course this model is more realistic than 
model B because of the smaller inner radius (accretion radius). 
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Fig. 5. Model C, high resolution (78 pc)and heating with cooling: the 
simulated gas density and temperature profiles are the same as in pure 
cooling simulation. The big dots represent the observations 



4.4. Model D 

In the model D, we set cutoff radius (ro - 44 kpc) at a larger 
value than in model B in order to get a smoother heating profile, 
see table 1 . The result is much better than model B but there is 
still a problem to fit the observations, see figure 6. The tempera- 
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Fig. 6. Model D, coohng with heating: the simulated gas density and Fig. 8. Model A, pure cooling: the cooling time scale is very short in the 
temperature profiles as functions of radius. The big dots represent the center, even there is no much cooling. 
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Fig. 9. Model B: the cooling time scale is higher than pure cooling sim- 
Fig. 7. Model D (r„ = 44 kpc): the plot of the ratio between the heating w t »i, » i - ,■ i ■ u- u ,u u j 

= , ,. ._,,„,.. . . . .. ulation. In the center, the coohng time scale is higher than observed 



rate and cooling rateCTYv / Ey). With increasing the the cutoff radius ro, 
the heating profile becomes more smoother. 



which is against the heating of cluster. 



ture at the center is higher than observed. With the large cutoff 
radius r,,, the heating profile can balance the cooling but it be- 
comes an artificial model which is not physical. In other words, 
the gradient of the heating profile becomes near to the gradient 
of the X-ray emissivity, ensuring the balance between the heat- 
ing and cooling. 

4.5. The effect of ttie tieating on cooling time scale 

In some cooling flow clusters, the observed cooling time scale 
is very short, in the range of 10** yr to 10'' yr, but the cooling 
gas below 1-2 keV is not observed. This result is interpreted as 
indication that there must be a heating mechanism to stop the 
cooling. In figure 8, The cooling time scale is very short in the 
center, even there is no catastrophic cooling or there is no much 
cooling. Figure 9 show that the heating increases the cooling 
time scale. If the cluster was heated, then we should find that the 
cooling time scale must be higher than observed. 



5. Summary and conclusion 

We have carried out time dependent simulations with cooling 
and heating. The general result is that the best fit to the observa- 
tions is the model without heating (model A). The effervescent 
heating is a function of the ICM pressure gradient but the 
cooling is proportional to the gas density square and square root 
of the gas temperature. Furthermore, the conclusions are : 

1- From our simulations, the cooling can not be a steady 
flow as assumed by standard cooling flow model. The gas 
density must increase with the time; i.e the ICM must be 
compressed under the force of the inflowing gas. 

2- The effervescent heating model profile (without the 
smoothness function) is more steeper than the radiative cooling 
profile, which makes it is very difficult to balance the cooling. 

3- The function 1 - e^''"^'' is not a cutoff function but it is 
a smoothing function in order to control the steepness of the 
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heating function letting it balance the radiative cooling; i.e. it is 
an artificial model. 

4- At inner radius close to the realistic value of the accre- 
tion radius (model D), the accretion mass rate due to flowing 
gas is not enough to fuel the black hole in the center of cluster. 
It's not reasonable to set the accretion radius of the black hole 
at large radius, (for example 1-1.5 kpc, as used for the three 
dimension simulations). 

5- In some cooling flow clusters, the observed cooling 
time scale is very short, in the range 10** yr to 10^ yr, but 
the cooling gas below 1-2 keV is not observed. This result is 
interpreted as indication that there must be a heating mechanism 
to stop the cooling. As in figure 8, the cooling time scale inside 
a radius of 10 kpc is shorter than 1 Gy but there is no a cool 
gas. Moreover, from our simulations (model B), we found that 
the heating process increases the cooling time scale. If cooling 
time scale is a good approximation for the actual cooling time 
and the cluster is heated, then we must find that the cooUng time 
scale is larger than observed. 
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